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Abstract: In this paper, we investigate thick branes with a non-minimally coupled back- 
ground scalar field, whose solution is a single-kink or a double-kink. The efi'ects of the 
non-minimal coupling constant on the structure of the thick branes and the localization 
of gravity, fermions and scalars are discussed. It is shown that each brane will split into two 
sub-branes as increasing the non- minimal coupling constant ^. By investigating the tensor 
perturbation equations of gravity and the general covariant Dirac equation of fermions, we 
find that both the gravity zero mode and left-chiral fermion zero mode are localized at 
the center of the single-kink branes and localized between the two sub-branes generated 
by the double-kink, which indicates that the constant ^ does not effect the localization of 
these zero modes. However, the zero mode of scalars is localized on each sub-brane (for 
both single-kink and double- kink branes) when ^ is larger than its critical value ^o- The 
effects of the non-minimal coupling constant ^ on the resonances of gravity and fermions 
with finite lifetime on the branes are also discussed. 
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1. Introduction 



The idea that our observed four-dimensional universe might be a hypersurface, embedded 
in a higher dimensional space-time (the bulk), provides new insights solving the gauge 
hierarchy and cosmological constant problems ||l|, |2|, ^, ^, ^, |7|, §]. In the Randall- 
Sundrum (RS) braneworld model, the zero mode of gravity is localized on the brane, which 
reproduces the standard Newtonian gravity on the brane |Q] . But in this model, the brane 
is very ideal because its thickness is neglected. In the most fundamental theory, there seems 
to exist a minimum scale of length, thus the thickness of a brane should be considered in 
more realistic field models. For this reason, more natural thick brane scenarios have been 
investigated g 0, |lT], |l2|, |l|, |l|, ||, |l6|, l^, ll 



19, pOL m\, 22, 23, 24, 25, 26, 27, 28, 29, 30 



3l|, H, H, H, H, H, |3> Hi' m ii- For some comprehensive reviews about thick branes, 



see Refs. |4l|, |2|, ||, |4 



In braneworld theory, the localization of gravity and various bulk matter fields is a very 
important issue. In order to recover the effective four-dimensional gravity, the gravity zero 
mode should be localized on branes. In the other hand, various bulk matter fields should 
be localized on branes by a natural mechanism, for the purpose to build up the standard 
model. Generally, the massless scalar fields |47| and graviton |5[ could be trapped on 
branes of different types. Spin-1 Abelian vector fields can be localized on the RS brane 



in some higher dimensional cases |48] or on the thick dS branes and Weyl thick branes 
1 54]. It is important to study the localization problem of the spin-1/2 fermions. Without 
introducing the scalar-fermion coupling, fermions can not be localized on branes in five and 
six dimensions f^, g ||, |0|, |l], ||, ||, ||, ||, |6|, ^ g ||, |0|, |l], ||, ||, ||, ^, ||, 



68, 69, 70 1 . In some cases, there may exist a single bound state and a continuous gapless 



spectrum of massive fermion Kaluza-Klein (KK) modes [^, 54]. In some other cases, one 
can obtain finite discrete KK modes (mass gap) and a continuous gapless spectrum starting 
at a positive [26, p^, |63[. 



Generally, including nonlinear terms of the various curvature tensors (Riemann, Ricci, 
Weyl) and non-minimally coupled terms in the effective action of gravity is a very com- 
mon trend from quantum field theory and cosmology. These theories cover f{R) modified 
gravity, the Gauss-Bonnet gravity, scalar-tensor gravity, and so on. Recently, braneworld 
models with a non-minimally coupled bulk-scalar field, via an interaction term — 
where (j) is a background scalar field and ^ is a non- minimal coupling constant, have been 
studied. Static solutions of these models have been examined in Ref. [71, 72, 73, 74 [. 
The Newton's law for brane models with a non-minimally coupled bulk scalar field was 
investigated in [|7H|. The conditions for localization of gravity for thick braneworlds with 



non-minimally coupled term and Gauss-Bonnet term were discussed in [76[. 

In this paper, we investigate thick branes with a non-minimally coupled background 
scalar, which is set as a single-kink and a double-kink. The structure of the thick branes is 
discussed, and the character of the branes is related to the non- minimal coupling constant 
^. For the single-kink scalar ^ becomes larger, the brane will split into two sub- 

branes. For the double-kink scalar case, there are two sub-branes located at both sides of 
z = 0, and as ^ becomes larger, the distance of the two sub-branes will increase. Further, 
the effects of the non-minimal coupling constant on the localization of gravity and various 
matter fields are investigated, and we find that the results especially for gravity, scalars 
and fermions are very nontrivial. For gravity and fermions, we find that the zero modes 
are localized on the center of branes generated by single-kinks or between the two sub- 
branes by double-kinks. For scalars, with the increase of ^, the zero mode is localized on 
the center of the two sub-branes or between them first and then on each sub-brane. The 
effects of the non-minimal coupling constant ^ on the resonances of gravity and fermions, 
which correspond to gravitons and fermions with a finite lifetime on the branes, are also 
discussed. 

The organization of this paper is as follows: In Sec. ^ the model of thick branes with 
non-minimally coupled bulk-scalar field in five-dimensional space-time is described. Then, 
in Sec. ^, we investigate the structure of the thick branes and localization of gravity on the 
branes. In Sec. ^, the localization of various bulk matter fields is investigated. Finally, the 
conclusion is given in Sec. ^. 



2. The model of thick branes with a non-minimally coupled bulk-scalar 
field 

We start with the following 5-dimensional action of thick branes, which are generated by 
a real non-minimally coupled scalar field (p, 



S = / d^x^ 



F{<l))R-]^g^''dM<t^dN<t^-V{<t^) 



(2.1) 
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where R is the five-dimensional scalar curvature and -F(</>) is chosen as 

F{ct>) = ^{l-4ict>''), (2.2) 

5 

= SvrGs with G5 the 5-dimensionaI Newton constant and is a dimensionless coupling 
constant. The 5-dimensional cosmological constant has been included in the scalar potential 

The Einstein equations corresponding to the action ( p.!]) are expressed as follows 

Rmn — i^QmnR = Tmn (2.3) 

with Tmn the energy- momentum tensor for the scalar field: 

'1 



Tmn = dM4>dN4> - 9mn 



2^ jp(l>dQct> + v{4>) 



+ 2VmV^F(0) - 2gMNUF{(t)) 



+(1 - 2F(0)) {rmn - \9mnR^ , (2.4) 

where □ is the five-dimensional D'Alembert operator. We have set = 1. 
The equation of motion for the scalar field reads 

, dF(d)) dVU) , , 



The above equation and the Einstein equations ( |2.3D are not independent ||7l[| . 

The line-element for the background spacetime describing a thick flat brane is assumed 

as 

ds'^ = gMNdx^dx^ = e'^^^y^T^f.^dx'^dx'' + dy^, (2.6) 

where e^"^^^^ is the warp factor, r/^j^ = diag(— 1, +1, +1, +1) is the Minkowski metric, and 
y stands for the extra coordinate. We suppose that the scalar field is considered to be 
a function of y only, i.e., cp = (j){y)- From Eqs ( |2.3D -(^^), we can obtain the following 
equations: 

3(1 - C<I>^){A" + 2A'^) + (^ - 2e)</.'2 + y(0) - 2^00" - 6^A'#' = 0, (2.7a) 

6(1 - ^cP^)A'^ - ^(A" + V{<P) - 8^^'#' = 0, (2.7b) 

- ct>" - 4A'cl)' - ^{8A" + 20^'2)</. + = 0, (2.7c) 

d(p 

where the prime stands for the derivative with respect to the extra coordinate y. A static 
analytical solution, for a narrow range of the coupling constant values < ^ < g, has been 



studied in Ref. |73|. In this paper, we will investigate the effect of the coupling constant 
^ on the brane structure and the localization of gravity and various spin fields in some 
general situations. 
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3. The structure of the thick branes and locahzation of gravity on the 
branes 



In this section we will investigate the inner structure of the brane with the non-minimally 
coupled background scalar field as well as the localization of the gravity on them, by the 
linearized equations for the metric fluctuations. From the Einstein equations ( |2.3| ), we can 
obtain the alternative form of the Einstein equations: 



2F{<P)Rmn = Tmn — -^QmnT, 



(3.1) 



where 



Tmn = dM4'dN(t> — Qmn 



\gPQdp<t>dQcl> + V{<l,) 



2eVAfV^<^ + 2igMNn^ (3.2) 



with T = q^'^^Tmn- From the above equations, we will obtain the linearized equations for 
the metric fluctuations Hmn- 

Under the axial gauge conditions h^M = 0, the total metric can be written in the form 



The Ricci tensor can be computed from the metric ( |3.3| 



RmN — -^MiV + ^MTV + 



where the zero order terms are 



0, 



(3.3) 



(3.4) 



(3.5) 



and the first order terms are 

rW = -e'^il-d'y + 2A'dy + A" + 4A'2) V - ^n^'^ V " \v,ue'^A'dy{r,'^^Kp), 



R, 



R, 



(1) 

55 

(1) 
5fi 



{dl + 2A'dy)r]'"'h^u, 



And we can also obtain 



Tmn - ^Ifiv + ^Mjv + 



with 



'f(o) 

fjLU 



n{0) 
'55 



Q—S'A 



+ 2^(b" + ^</.'2 + y(0) 



1 



<^'2 _ y(^) _ 8 



p{0) 



0, 



(3.6) 
(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 
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and 



dF ,, (fF 



1 



+ 2'^(t)" + + V{(t)) 



(3.12) 
(3.13) 



where h = rj^'^h^y. From above equations and considering the fluctuations h^y satisfy the 
transverse d^h^y = and traceless h = rj^^h^y = conditions, we can obtain the foUowing 
equations 



dl + Q'{y)dy + e-2^n(4)l hl^ix, v) = 0, 



(3.14) 



where we have set 



Q{y) = 4.A{y) + \nF{<l){y)). 



(3.15) 



Here, F{(j){y)) should satisfy F{<j){y)) > 0. By using the decomposition hj^i^{x,y) 
e'^P^hfj_y{y), the above equation can be reexpressed as 



[d^ + Q'iy)dy + m\-^^] V(y) = 



(3.16) 



with p'^p^ = —m? the four-dimensional mass of the gravitons. 

For the sake of convenience of obtaining the mass-independent localization potential 
for gravitons, we will follow Ref. [^] and change the metric given in (|2.6| ) to the conformally 
flat form 



^ Q^M-^){^^^dx^'dx'' + dz"^) 
by performing the coordinate transformation 

dz = e-^^y^y. 

Hence Eq. (|3.14| ) can be rewritten as 

d^,+Q'{z)d,+m^] V(z) = 0, 



(3.17) 



(3.18) 



(3.19) 



where Q{z) = 3A{z) + lnF{(p{z)). By using the transformation h^i,{z) = e ^/'^Lp^y{z)^ Eq. 
( 3.191) can be expressed as the Schrodinger equation 

[-dl + V^uiz)] ^{z) = mV(^), (3.20) 

where ip^y{z) has been replaced by (p{z) for simplicity and the localization potential is read 



as 



^qm(^) = \q"{z) + \q'\zI 



(3.21) 
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where the prime denotes the derivative with respect to z. Eq. ( |3.2C1| ) can be rewritten 
alternatively as 



9. 



(3.22) 



so it will have a complete system of eigenstates with non-negative eigenvalues, i.e., w? > 0. 
By setting m = 0, the massless wave function (the zero mode) can be obtained 



99o(^) ocei^(^)F5((/,(z)). 



If the zero mode satisfies the normalization condition 

f + 00 



(pQ{z)dz 



33^(")F(</.(z))dz 



< oo, 



(3.23) 



(3.24) 



the zero mode is localized on the brane. 

From above discussion, it is easy to see that static solutions for the thick branes are 



convenient for calculation in the z coordinate with metric ( 3.17 ). Hence, we will study the 
static solutions in the z coordinate. The equations of the motion for the background scalar 
field (j){z) and the warp factor A{z) in the z coordinate can be written as 

3(1 - e<^2)(^'2 ^ ^//) ^ 1^.2 ^ e2^y((/>) - 2CcP'^ - 2^cPcP" - ^i4>4)" = 0, (3.25a) 

6(1 - C</)2)^'2 - i,/,'2 + e2^y(0) - 8e0^'0' = 0, (3.25b) 

^" + + 4e</.(3^'2 + 2^1") - e^^^^ = 0. (3.25c) 

dcp 

As we have already mentioned, above equations are not independent. Hence, we choose to 
solve Eq. ( p. 25c ) and the following equation, which can be obtained by ( |3.25a ) minusing 



(|3]25bD, 

3(1 - i4>'^){A" - A'"^) + Ai4)A'4)' + (1 - 2e)</)'2 - 2i4)(t)" = 0. 
In this paper, the background scalar field is considered as a kink solution 

(j){z) = 0otanh'''(^^;), 



(3.26) 



(3.27) 



where and b are positive real parameters, and the parameter k is an positive odd number. 
Here we just only consider two cases: the k = 1 single-kink solution and the A; = 3 double- 
kink solution. The background scalar has the following behavior at z = and z — >• ±00: 



m = 0, 



±00 



(3.28) 



Because -F((/>) > and (j) < 0o, the coupling constant ^ < It is very difficult to 
solve Eq. (|3.26D by analytical method, so we will numerically solve it. First, we will 
analyze the asymptotic behavior of the warp factor at 2 — ?■ ±00. When z — >• ±00, we have 
for the scalar field that (j) — >• ^(f^o, cj)' ^ and (p" — 0. Hence, Eq. ( |3.26D reduces to 
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A"[z) — A''^{z) = as z — 7- ±00, from which the asymptotical form of the warp factor can 
be given as e^*-^-* = 1 ~ m, which is the same as the RS braneworld. Next, we want 



to obtain the approximate solution of Eq. (3.26) in the vicinity oi z = 0. We expand the 
background scalar field (j){z) at 2: = 0: 

^{z) = Mbz)' - l(l^ok{bz)>'+^ + 0[hzf^\ (3.29) 

For = 1, by using i?i>(z) = 4>o{bz) — ^(f)o{bz)^ + 0{bz)^, the behavior of the warp factor 
in the vicinity of z = 0, from Eq. ( p. 26 ), can be expressed as 



A{z) = ^4>1{2C - l){bzf + 0{bzf. (3.30) 

We find that the behavior of the warp factor in the vicinity of z = is related closely to 
a critical coupling constant ^0 = ^- When small ^ with ^ < ^O) the maximum of the warp 
factor e^'^(^) is at z = 0, and this will correspond to one brane located at z = 0. When ^ is 
closer to we will see that the thick brane will split into two sub-branes. When > (,o, 
the maxima of the warp factor are at both sides of z = and this will result in the increase 
of the distance of the two sub-branes. 

For k = 3, the scalar 4>{z) can be expressed as (j){z) = 4>o{bz)'^ + 0{bz)^ at z = 0. And 
the warp factor can be described as 

A{z) = ^<Plm - mzf + o{bzy. (3.31) 

Here the critical coupling constant is .^o = The behavior of the warp factor is the same 
as the case k = 1. However, there are always two sub-branes for any ^ because the scalar 
is a double-kink. 

Eq. ( |3.26| ) can be solved numerically with the following initial conditions: 

A(0) = ^'(0) = 0. (3.32) 

The shapes of the warp factor e^"^ and the energy density Too(z) are shown in Fig. || and 
Fig. 1^ for k = 1 and /c = 3, respectively. For k = 1, when ^ is very small, there is only one 
brane located at z = 0, and when increasing ^, the brane gradually splits two sub-branes 
located at both sides of z = 0. The distance of the two sub-branes increases with ^. For 
A; = 3, there are two sub-branes located at both sides of z = even for small ^, and the 
distance of the two sub-branes also increases with ^. 

From these solutions, the potential of gravity KK modes Vqm (|3.21| ) and the gravity 
zero mode ( 3.23| ) are also can be obtained, which are shown in Figs. ^ and ^ 



By considering the asymptotic behavior of the warp factor and the background scalar, 
the asymptotic behavior of the potential as z — t- iLcxd is 

Vqm{z ±oo) ^ 0, (3.33) 

which is similar to a volcano-type potential. The asymptotic behavior of the gravity zero 
mode is 



^o{z ^ ±oo) oc ^(1 - ^</.g)i ' . (3.34) 



3 
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(a) = 1 (b) = 1 

Figure 1: The shapes of the warp factor e^^ and the energy density Too (2). The parameters are 
set to ^ = 0.1 for the thiclc red line, ^ — 0.5 for the dashed green line, and ^ = 0.9 for the thin blue 
line. The other parameters are set as c/jq = 1, b = 1, and fc = 1. 



^2 At) 

7"oo(z) 




(a) = 3 (b) = 3 



Figure 2: The shapes of the warp factor e^ and the energy density Too(z). The parameters are 
set to ^ = 0.1 for the thick red line, ^ = 0.3 for the dashed green line, and ^ = 0.9 for the thin blue 
line. The other parameters are set as c/jq ~ 1, b ~ 1, and fc = 3. 



and the normalization condition ( p. 24 ) is equivalent to the following condition 



so it is clear that the normalization condition is satisfied. For k = 1, when ^ is smaller than 
the critical coupling constant the zero mode is localized on the single brane, and when 
^ is close to (,0 or ^ > S^Q, the zero mode is localized at the center of the two sub-branes. 
For k = 3, the zero mode is localized between the two sub-branes. The massive gravity 
KK modes are the continuous spectrum wave functions with > 0, and asymptotically 
turn into plane waves. 
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(a) = 1 (b) A: = 1 

Figure 3: The shapes of the potential of gravity KK modes Vqm(^) and the gravity zero mode 
ipoiz). The parameters are set to ^ = 0.1 for the thick red hne, ^ = 0.5 for the dashed green hne, 
and ^ — 0.9 for the thin bhie hne. The other parameters are set as (/)o = Ij & = 1, and k = 1. 




(a) it = 3 (b) A: = 3 

Figure 4: The shapes of the potential of gravity KK modes Vqm{z) and the gravity zero mode 
(Pq{z). The parameters are set to ^ = 0.1 for the thick red line, ^ = 0.3 for the dashed green line, 
and ^ = 0.9 for the thin blue line. The other parameters are set as (/iq = 1, & = 1, and fc = 3. 



Generally, the volcano-type potential implies that there may exist resonant states, 
which are quasi- localized on the branes. So we will study whether the resonant states exist 
in this system. From Figs. ^ and ^, it is seen that with the coupling constant ^ becoming 
larger, the height of the potential Vqm gets higher, with which the resonant states could 
be obtained more easily. When the coupling constant ^ = 0.9 and other parameters are set 
as 00 = 1 and 6 = 1, there is only one resonant state for the double- kink background scalar 
field. However, we find that when the coupling constant ^ = 0.99 and other parameters are 
set as 00 = 1 and 6 = 1, there is a resonant KK mode for the single-kink background scalar 
field, and two resonant KK modes for the double-kink background scalar field. Because 
resonant states are oscillating when far away from the brane along the extra dimension. 
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they can not be normalized. As Ref. |59], we use the function 



(3.36) 



to describe the relative probability for finding the resonances on a thick brane. Here 2z\) is 
about the width of the thick brane and z^^x is set to z^^x = lOz^. It is clear that for KK 
modes with rr? » ^qm^ (^qm^ maximum value of Vqm), which are approximately 

taken as plane waves, the value of PQ,{rr?') will trend to 

The profiles of Pq corresponding to the coupling constant ^ = 0.9 for single-kink scalar 
field (fc = 1) and double-kink background scalar field (k = 3) are shown in Fig. |5| for the 
thick brane, and Pq corresponding to ^ = 0.99 for k = 1 and A; = 3 are shown in Fig. In 
these figures, each peak corresponds to a resonant state. When ^ = 0.9, it is seen that only 
for the double-kink background scalar field there is a resonant KK mode. When ^ = 0.99, 
there are one and two resonant KK modes for single-kink and double-kink background 
scalars, respectively. Hence, we come to the conclusion that if the coupling constant ^ 



becomes larger, the number of the resonant gravity KK modes will increase. As Ref. jl5|, 
we estimate the lifetime r of a resonant state as r ~ F"^ with F = 6m being the full width 
at half maximum of the peak. Then the mass m, width F and lifetime r of the resonant 
KK modes are listed in Table ||, and the resonant KK modes are shown in Figs. |7|, |8| and 





(a) C = 0.9, k 



(b) C = 0.9, k = 3 



Figure 5: The profiles of Pq for massive even and odd gravity KK modes for the parameters 
^ = 0.9, </«o = 1, = 1, fc = 1 and fc = 3. 



4. Localization of the various matters on the thick branes 

In this section, we will investigate the effect of non-minimally coupled parameter ^ on the 
localization and resonances of various bulk matter fields on the thick branes, by presenting 
the potential of the corresponding Schrodinger equation for the KK modes of the various 
matter fields. Spin-0 scalars, spin-1 vectors and spin-1/2 fermions will be considered by 
means of the gravitational interaction. 
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(a) C = 0.99, k = 1 



(b) C = 0.99, k = 3 



Figure 6: The profiles of Pq for massive even and odd gravity KK modes for the parameters 
£, = 0.99, 00 = 1, 6 = 1, A: = 1 and fc = 3. 





k 


Height of Vqm 


n 


2 

m 


m 


r 


r 


0.9 


3 


y-- = 1.1322 


1 


0.685649 


0.82804 


0.123867 


8.07321 


0.99 


1 


V^^^ = 2.8051 


1 


1.634424 


1.27845 


0.041632 


24.02015 




3 


V^^^ = 2.7746 


1 


0.831835 


0.91205 


0.013373 


74.77949 








2 


2.433646 


1.56001 


0.104536 


9.56606 



Table 1: The mass, width, and lifetime of the resonances of gravity. The parameters are set to 
(po = 1 and 6=1. Here n is the order of the resonant states with corresponding from small to 
large. 



4.1 Scalar fields and vector fields 

We first consider the localization of scalar and vector fields on the thick branes obtained 
in the previous section, then turn to fermions in the next subsections. Let us start by 
considering the action of a massless real scalar coupled to gravity 



Sn 



(4.1) 



By using of the KK decomposition <I>(x, z) = (j)n{x)xn{z)e '^"^Z^, it is easy to derive the 
equations for the scalar KK modes: 



[-dl + Voiz)] Xn{z) = mlxniz), 
where the Schrodinger potential is given by 

Vo{z) = ^A"{z) + Ia'\z). 



(4.2) 



(4.3) 
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1.0 



0.5 




10 



20 



(a) i = 0.9, A; = 3, n = 1 



Figure 7: The wave function of the gravity resonant KK mode for the double-kink background 
scalar. The parameters are set to ^ = 0.9, (/jq = 1, 6 = 1 and fc = 3. 



Figure 8: The wave function of the gravity resonant KK mode for the single-kink background 
scalar. The parameters are set to ^ = 0.99, (f){) — \,h — \ and fc = 1. 

Here rrin are the masses of the scalar KK modes, they are also the masses of the four- 
dimensional scalars (/)n{x). It is clear that the potential Vo{z) defined in ( [4. 3D is a four- 
dimensional mass-independent potential. Note that those scalar KK modes localized on 
the branes should satisfies the following orthonormality conditions: 



For the thick branes considered previous section, the warp factor can not be written 
as an analytical form, however, the potential Vq can be solved numerically, which is shown 
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(a) 5 = 0.99, A: = 1, n = 1 




(4.4) 



in Figs. ^ and 11. 



The scalar zero mode can be solved from (|4.2D by setting uiq = 0: 



(4.5) 
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(a) ^ = 0.99, A; = 3, n = 1 



(b) C = 0.99, fc = 3, n = 2 



Figure 9: The wave functions of the gravity resonant KK modes for the double-kink background 
scalar. The parameters are ^ = 0.99, = 1, 6 = 1 and A: = 3. 



This scalar zero mode (hence the four-dimensional massless scalar) is localized on the 
branes because the boundaries of the five-dimensional spacetime along extra dimension are 
AdS. This can also be confirmed by the following simple calculation: 



Xo{z 



1 



±oo) 

3 



1 



dz 



1 

263 



< oo. 



(4.6) 
(4.7) 



In fact, this is a well known conclusion in braneworld model. Here, our point is another 
question: where the four-dimensional massless scalar locates when the thick branes split? 
As far as we know, when a thick brane splits into two sub-branes in minimally coupled 
theories, the four-dimensional massless scalar is located between the two sub-branes. Here, 
for the case of non-minimally coupled theory considered in this papaer, we can see from 
Figs. 10 and 11, as well as Eqs. (4^), ( ^3.30 ) and ( 3.31 ) that, the four-dimensional massless 
scalar is localized at the center of the two sub-branes generated by single-kink or between 
them by double-kink as ^ < just as the usual situation, while it is localized on each 
sub- brane when the coupling parameter ^ is greater than its critical value ^o- 

All other scalar KK modes are continuum and massive, for which we do not find 
resonant KK modes in the spectrum. This is different the case of gravity. 

For Abelian spin-1 vectors described by the following action 



Si 



— MN RS TP 77 

-g g g FmrFns, 



(4.8) 



the zero mode can not localized on the branes with AdS boundaries. This problem can be 



solve in some higher dimensional models |48] or in thick dS branes and Weyl thick branes 



models |54]. We do not yet find vector resonant KK modes in the spectrum. 
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(a) = 1 (b) A: = 1 

Figure 10: The shapes for the potential of the scalar KK modes Vq{z) and the scalar zero mode 
Xo{z)- The parameters are set to ^ = 0.1 for the thick red line, ^ = 0.5 for the dashed green line, 
and ^ = 0.9 for the thin blue line. The other parameters are set as = 1, ^ = 1, and k — 1. 




(a) it = 3 (b) it = 3 

Figure 11: The shapes for the potential of the scalar KK modes Va{z) and the scalar zero mode 
Xo{z)- The parameters are set to ^ = 0.1 for the thick red line, = 0.3 for the dashed green line, 
and ^ = 0.9 for the thin blue line. The other parameters are set as (j>o ^ 1, b — 1, and k ^ 3. 



4.2 Spin-1/2 fermion fields 

Next, we will study the localization and resonances of fermions on the thick branes. In 
five-dimensional space-time, fermions are four-component spinors and their Dirac structure 
can be described by the curved space gamma matrices T^'^ = e~^{'y^,j^), where and 
7^ are the usual flat gamma matrices in the four-dimensional Dirac representation. The 
Dirac action of a massless spin-1/2 fermion coupled to the background scalar (p ( ^.27 ) is 



Si/2 = j dJ'x^ {^>T^'D^^ - r?^G((/.)*) , (4.9) 

where r/ is a coupling constant, and D^"^ = {Bm+^m)'^ with the non-vanishing components 
of the spin connection ojm for the background metric ( |3.17 ) given by = ^^'7^75. Then 
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the equation of motion reads 



[rd^ + 7' {d, + 2 A') - 77 e^G{ct>)\ = 0, (4.10) 

where ^^d^ is Dirac operator on the brane. 

Now we would hke to investigate the effect of the couphng parameter ^ on the lo- 
cahzation and resonances of the Dirac spinor on the branes by studying the above five- 
dimensional Dirac equation. Because of the Dirac structure of the fifth gamma matrix 
7^, we expect that the left- and right-chiral projections of the four-dimensional part have 
different behaviors. From the equation of motion ( [4.1[)| ), we will search for the solutions of 
the general chiral decomposition 

^ = 2^ [V'L,n(x)L„(z) + Vi?,n(x)i?„(z)] e'^^, (4.11) 
n 

where ipL = V' s-^id = ^-y-V' ^-re the left- and right-chiral components of a four- 
dimensional Dirac field respectively. By demanding V'L,/? satisfy the four-dimensional 
massive Dirac equations '^^{d^ + u}^)i/jl,r = fn'4'R,L-, we obtain the following Schrodinger- 
like equations for the KK modes of the left- and right-chiral fermions: 

{-dl + Vl{z))L^ = n?K, (4.12a) 
{-dl + VR{z))Rn = m^Rn, (4.12b) 

where the effective potentials are 

Vl{z) = {r]e^G{^))^ - d^ve^Gi^)), (4.13a) 

Vr{z) = VL{z)\r,-,-r,. (4.13b) 

For the purpose of getting the standard four-dimensional action for a massless fermion 
and an infinite sum of the massive fermions: 

5i = y d^xV^ ^ [T^'idM + com) - r]G{^)] ^ 

= / d^^"/-^ Ml'^id/i + (^t^) - mn]i^n, (4.14) 
we need the following orthonormality conditions for L„ and Rn- 

LrnLndz = 6mm (4-15) 
RmRndz = 6mn, (4.16) 

LmRndz = 0. (4.17) 



oo 



oo 
oo 



-oo 
oo 



From Eqs. (p^) and (|1|), it is clear that, in order to localize the left- and right- 
chiral fermions, some kind of scalar-fermion coupling must be introduced. This situation 
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is similar to the one in Refs. |47, 48, 54, 20, 51, 59, in which the authors introduced 



the scalar-fermion coupUng term m^F((j))^ for the locahzation of the fermion fields on 
a brane. Moreover, if we demand that Vl{z) and Vr{z) are Z2-even with respect to the 
extra dimension coordinate z, G{(j)) should be an odd function of the kink (j){z). In this 
paper, we choose the simplest Yukawa coupling: G{(j)) = (f). So the potentials for left- and 
right-chiral fermion KK modes can be expression as 



VLiz) 



2 2A(z)i2/ \ 

f] e ^ '(p [z) 



(4.18) 
(4.19) 



First, by considering the initial conditions of the warp factor and the background scalar at 
z = 0, the values of the potentials can be obtained: 



Vl{0) 
Vr{0) 



-r](j)'{z) 
+#'(z). 



z=0 ' 
z=0 ' 



So for k = 1, 



and for k 



Vl{0) = -#o6, Vjj(O) = +7j<Pob, 



VlA^) = 0. 

Second, when far away from the brane z — )• ±00, 

Vl,r{z ^ ±00) ^ 0. 



(4.20a) 
(4.20b) 



(4.21) 



(4.22) 



(4.23) 



Here the potentials can be solved by numerical method, which are shown in Figs. H and 



13. 



The left- and right-chiral fermion zero modes are solved as 



Lo{z) oc exp ( -7? / dz'e^^^'Uiz') 



Ro{z) oc io(^)U- 



-v 



From the asymptotic behavior of the warp factor as 2: — )• 
left- and right-chiral fermion zero mode can be analyzed: 



Lo{z — ±00) 
Ro{z — )• ±00) 



\z\ b 
I I H^o. 

\Z\ b . 



(4.24a) 
(4.24b) 

±00, the asymptotic behavior of 

(4.25a) 
(4.25b) 



So for the positive coupling constant r/, only the left-chiral zero mode tends to zero when 
far away from the branes, which may be localized on the brane. We need further check 
whether the normalization condition ( 4.15| ) is satisfied for the left-chiral zero mode, i.e., 



dzL^{z) < 00. 



(4.26) 
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Since the values of the zero modes are finite at finite z, the above normahzation condition 
is equivalent to the following condition 



2^700 

dz\z\ — — < oo. (4.27) 

Only when r] > tjo = (note that 00 and b are positive real parameters), the above 
integral is convergent, which means that the left-chiral zero mode can be localized on the 
branes under this condition. From the shapes of the zero mode of left-chiral fermions 



(Fig. 14), it can be seen that the zero mode is localized on the center of the two sub- 
branes when the thick brane splits (single- kink case) or between them (double- kink case). 
However, we will show in the following that the non-minimally coupling parameter ^ would 
eff'ect the resonance spectrum of both the left- and right-chiral fermions. 




(a) it = 1 (b) it = 1 

Figure 12: The shapes of the potentials for the left- and right-chiral fermions coupled with single- 
kink scalar (fc = 1). The parameters are set to ^ = 0.1 for the thick red line, ^ = 0.5 for the dashed 
green line, and ^ = 0.9 for the thin blue line. The other parameters are set to 0o = 1, r; = 5, and 
6=1. 



All the massive KK modes of the left- and right-chiral fermions are the continuum 
modes, and can not be localized on the brane. Next, we are going to investigate the quasi- 
localization of the left- and right-chiral fermion KK modes, which are called resonances. 
As Sec. ^ the relative probability for finding the resonances on a brane can be given as 

r" \L,R{z)\^dz 

In order to investigate the effect of the gravity-scalar coupling constant on fermion reso- 
nances, we will investigate the resonant KK modes with different values of ^. 

For the single- kink background scalar field {k = 1), the mass, width and lifetime of the 
resonant KK modes with different values of ^ are listed in Table |2|. As an example, we plot 
the shapes of -Pl,r(?7i^) corresponding to ^ = 0.9 and A: = 1 in Fig. |l5[ From Table |2| and 
Fig. f5|, it can be seen that the mass and lifetime of left- and right-chiral fermion resonances 
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(a) A; = 3 (b) = 3 

Figure 13: The shapes of the potentials for the left- and right-chiral fermions coupled with double- 
kink scalar (fc = 3). The parameters are set to ^ = 0.1 for the thick red line, f = 0.3 for the dashed 
green line, and ^ = 0.9 for the thin blue line. The other parameters are set to = Ij = 5, and 

b=l. 



^0 La 




-2-1 12 -2-1 12 



(a) = 1 (b) A; = 3 

Figure 14: The shapes of the left-chiral fermion KK zero mode. The parameters are set to ^ = 0.1 
for the thick red line, ^ = 0.5 for the dashed green line, and ^ = 0.9 for the thin blue line for the left 
figure (fc = 1), and the parameters are set to ^ = 0.1 for the thick red line, ^ = 0.3 for the dashed 
green line, and ^ = 0.9 for the thin blue line for the right figure (fc = 3). The other parameters are 
set to 4>o = 1, rj — 5, and b — I. 



are almost the same, so the formation of the four-dimensional massive Dirac fermions can 
be realized |^8[. The shapes of the resonances Ln and Rn are shown in Figs. 16 and ^ for 
^ = 0.9 and k = 1. For the left-chiral KK resonances, it can be seen that the first resonant 
KK mode is an odd-parity wave function, and the second one has even parity. However, for 
the right-chiral KK resonances, the first resonant KK mode has even parity and the second 
one has odd parity. This is held for any n-th fermion resonances, namely, the parities of 
the n-th left- and right-chiral resonances are opposite. In fact, this conclusion is originated 
from relationship between the two potentials Vl{z) and Vii{z), or equivalently, from the 
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coupled equations of the left- and right-chiral fermions, which are not given in this paper 
but one can refer [^]. Further, four-dimensional massive Dirac fermions can be obtained, 
which are consisted of the pairs of coupled left- and right-chiral KK modes with different 
parities. 

From Table g it is shown that the lifetimes of the first resonant KK modes satisfy 
Tn=i{i = 0.5) > Tn=i{i = 0.1) > r„=i(^ = 0.9). From Fig. p], it is seen that the width 
of the barrier of the potentials in the vicinity of Vl,r = ml for ^ = 0.5 is close to that for 
^ = 0.1 and larger than that for ^ = 0.9, so Tn=i{C = 0.5) and Tn=i{^ = 0.1) are larger than 
Tn=i{C = 0.9). Moreover, the widths of the barrier of the potentials for ^ = 0.5 and ^ = 0.1 
are the same, but the height of the barrier for = 0.5 is larger than that for ^ = 0.1, so 
T„=i(^ = 0.5) > Tn=i{(, = 0.1). For the second resonant KK modes, the situation is similar 



to the first ones. As Fig. 12 shown, the width of the barrier of the potentials in the vicinity 
of Vl,r = m-l for ^ = 0.5 is the largest, so t„=2(? = 0.5) is the largest one. Hence, we can 
come to the conclusion that the lifetimes of the resonant KK modes are decided by the 
width and the hight of the barrier of the potential and the width is more important than 
the hight. 





Chiral 


Height of Vl,r 


n 


2 

m 


m 


r 


r 


0.1 


Left 


ymax ^ U Q^Q 


1 


8.5206 


2.9190 


0.001222 


818.64177 








2 


13.6638 


3.6965 


0.111288 


8.98573 




Right 


ymax ^ U_QJ2 


1 


8.5206 


2.9190 


0.001216 


822.65478 








2 


13.6463 


3.6941 


0.111813 


8.94354 


0.5 


Left 


ymax ^ IQ QH 


1 


8.8961 


2.9826 


0.000725 


1380.15215 








2 


15.1313 


3.8899 


0.063637 


15.71416 




Right 


y^a^ = 16.651 


1 


8.8961 


2.9826 


0.000722 


1385.67372 








2 


15.1188 


3.8883 


0.064661 


15.46516 


0.9 


Left 


ymax ^ Ig ^H 


1 


9.2353 


3.0390 


0.003074 


325.2705 








2 


16.3330 


4.0414 


0.081694 


12.2408 




Right 


= 18.466 


1 


9.2353 


3.0390 


0.003120 


320.4962 








2 


16.3363 


4.0418 


0.084267 


11.8671 



Table 2: The mass, width, and lifetime of resonances for fermions. The parameters are set to 
A; = l,?7 = 5, 00 = 1 a-nd 6=1. Here n is the order of resonant states with corresponding from 
small to large. 

For double-kink background scalar field {k = 3), because the mass, width and lifetime 
of the right-chiral fermion resonant KK modes is the same as that of left-chiral fermion 
resonances, we only list the mass, width and lifetime for the left-chiral fermion resonant 
KK modes with different values of ^ are listed in Table ^ It is known that the number 



-19- 



1 J- 

0.8 - 

0.6 - 

0.4 - 

0.2 - 

- 



odd 
even 







10 





(a) ^ = 0.9, k = l, Left-chiral 



(b) ^ = 0.9, k = l, Right-chiral 



Figure 15: The profiles of Pl.r, f^or tlie left- and right-chiral fermion resonances for the brane with 
the parameters = 0.9, 77 = 5, 0o = 1 and 6=1. 
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(a) n = 1 



(b) n = 2 



Figure 16: The profiles of the wave function for the left-chiral fermion resonances L„ with the 
parameters ^ = 0.9, rj = 5, (f>o = 1 and 6 = 1. 

of the resonant KK modes will increase when the coupling constant ^ becomes larger. It 
is also shown that the lifetimes of the first resonant KK modes satisfy rn=i(^ = 0.1) > 
Tn=i{^ = 0.3) > T„=i(^ = 0.9), which is because that the widths of the barrier of the 
potentials in the vicinity of Vl,r = ml for ^ = 0.1 are the largest and that for ^ = 0.9 are 
the smallest from Fig. For the other resonant KK modes, the situation is similar to 
the case of A: = 1. 



5. Conclusions and discussions 

In this paper, we first review the model of the thick branes with a non-minimally coupled 
background scalar field and then investigate the structure of the branes. In our model, 
the background scalar field is set as single-kink (k = 1) and double-kink {k = 3) solutions 
respectively. The non-minimally coupling between the gravity and the background scalar 
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(a) n = 1 (b) n = 2 

Figure 17: The profiles of tiie wave function for tlie riglrt-cliiral fermion resonances i?„ with, the 
parameters ^ = 0.9, = 5, = 1 and 6=1. 
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2.92 X 10"^° 
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8.8214 
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1.45 X 10"^ 


6.89 X 10^ 
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13.3115 


3.6485 


0.018426 


54.2705 
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Left 


ymax ^ ^^-^^^^gQ 
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3.2171 


1.79363 


3.91 X 10-1° 


2.56 X 10^ 








2 


8.9598 


2.9933 


8.18 X 10-^ 


1.22 X 10^ 








3 


13.7770 


3.7117 


0.010634 


94.04994 


0.9 


Left 


ymax ^ 20.395 


1 


3.2735 


1.8093 


9.03 X 10-6 


1.11 X 10^ 








2 


9.3293 


3.05439 


0.00114 


878.85943 








3 


14.9016 


3.8603 


0.01752 


57.07182 








4 


19.8003 


4.44975 


0.128982 


7.75318 



Table 3: Tlie mass, width, and lifetime of resonances for left-chiral fermion. The parameters are 
set tofc = 3,?y = 5, 00 = 1 and h = \. Here n is the order of resonant states with corresponding 
vn? from small to large. 



field is introduced via a term ^^Rcj)"^, and we find that the behaviors of the warp factor 
and the branes are related to the non- minimal coupling constant ^. When the non- minimal 
coupling constant ^ is smaller than its critical value ^o, the maximum of the warp factor 
is at z = 0, and when ^ > the maxima of the warp factor are at both sides of z = 0. 
For single-kink (k = 1), the brane will split into two sub-branes with the increase of the 
non-minimal coupling constant ^, and the distance of the two sub-branes increases with ^. 
For double-kink (fe = 3), there are two sub-branes, and the distance of them also increases 
with ^. Further, we investigate the effects of the non-minimal coupling constant on the 
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localization of gravity and various bulk matter fields on the branes. 

Firstly, the localization of gravity is considered. It is found that, for the case of single- 
kink (A; = 1), the gravity zero mode is localized on the single brane with small ^, and it 
is localized on the center of the two sub-branes with large ^. For the case of double-kink 
{k = 3), the gravity zero mode is localized between the two sub-branes. All the massive 
modes are continuous spectrum wave functions and could not be localized on the branes. 
However, for larger coupling constant ^, there could exist gravity resonant states, and the 
number of the resonances increases with ^. 

For scalar field, the zero mode is localized on the center of the branes {k = 1 case) 
or between them {k = 3 case) when ^ < ^Oj while it is localized on each sub-brane when 
^ > ^0- For the vector field, the zero mode can not be localized on the branes. All the 
massive modes for scalars and vectors are continuous spectrum wave functions and can not 
localized on the branes. There is not exist the resonant state for both fields. 

For spin-1/2 fermion fields, in order to localized the left- and right-chiral fermions, we 
introduce the usual Yukawa coupling rj^cj)^. We find that for positive Yukawa coupling 
constant 77 larger than its critical value 770, the left-chiral fermion zero mode can be localized 
on the branes. For the case of = 1, the left-chiral fermion zero is localized on the single 
brane with small £^ and on the center of the two sub-branes with large ^. For the case of 
A; = 3, the left-chiral fermion zero is localized between the two sub-branes. The massive KK 
modes asymptotically turn into continuous plane waves when far away from the branes. 
It is interesting that the well of the potentials for left- and right-chiral fermions becomes 
deeper and deeper with the increase of the non-minimal coupling constant ^, which leads 
to a series of massive fermions with a finite lifetime on the branes. The spectra of left- 
and right-chiral fermion resonances are the same, which demonstrates that a Dirac fermion 
with a finite lifetime on the branes can be composed from the left- and right-chiral fermion 
resonant KK modes. 
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